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Abstract: 

We calculate the radiative corrections to the nonleptonic inclusive B decay mode b —>■ cud 
taking into account the charm quark mass. Compared to the massless case, corrections 
resulting from a nonvanishing c quark mass increase the nonleptonic rate by (4-8)%, de- 
pending on the renormalization point. As a by-product of our calculation, we obtain an 
analytic expression for the radiative correction to the semileptonic decay b — > utu taking 
into account the r lepton mass, and estimate the c quark mass effects on the nonleptonic 
decay mode b — > ccs. 
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1 Introduction 



The past few years have witnessed a continuous interest in the physics of heavy quarks, 
in particular B meson decays. In addition to their potential role as precision tests of the 
Standard Model, they also provide important information on the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements Vet, Vbu, and others. From a theoretical point of view, 
inclusive widths tend generally to be cleaner than exclusive ones. This is in particular true 
for the decays of heavy hadrons, which are essentially short distance processes since the 
amount of energy released is large compared with a typical confinement scale of ~ 1 GeV. 
Indeed it can be shown, cf. [|l], |^, that up to corrections of order 1/ml, ui}, being the 
b quark mass, inclusive decay rates of beautiful hadrons are essentially the same as for 
the underlying quark processes and thus can reliably be calculated in perturbation theory. 
Moreover, the corrections suppressed by powers of the heavy quark mass can be evaluated 
in a systematic way, using the operator product expansion (OPE) technique combined 
with the heavy quark expansion [|l], |^. At present, it is widely believed that inclusive 
decay widths — semileptonic, rare radiative or hadronic — may be reliably predicted to 
an accuracy limited mainly by the unavoidable shortcomings of perturbative calculations 
such as scale- and scheme-dependence. There is little controversy that these calculations 
rest on a firm theoretical foundation. 

A high accuracy thus requires the calculation of radiative corrections to inclusive widths, 
a task pioneered by the authors of 0, where the radiative correction to the semileptonic 
decay width of the D meson was obtained from earlier calculations |^ of the muon decay in 
QED. At present, 0{as) corrections are known for both the semileptonic and the hadronic 
decay width [H, H, in the approximation of massless final state quarks. 

Taking into account the masses of the final state particles has proved to be a difficult 
task. A significant effort in this direction has been made by Hokim and Pham 0, who 
calculated one-gluon corrections to inclusive weak decays taking full account of all the 
masses. Nir [Q has further found an analytical expression for the radiative correction to 
the semileptonic h cev decay, the numerical results being given in Unfortunately, the 
results of [0 are not sufficient for the calculation of inclusive widths with renormalization 
group improvement to two-loop accuracy, since in addition to the correction calculated 
by Hokim and Pham one must take into account terms proportional to In M^/m^ which 
are equally important and are also in general mass-dependent. Thus, in practice of data 
analysis quark mass corrections to nonleptonic decays are generally being omitted. On the 
other hand, it is known that the effect of non-zero masses on the radiative corrections can 
be very large, especially if they involve colour-attraction in the final state. An example 
was provided recently within the heavy quark effective theory: radiative corrections to 
correlation functions, containing one heavy and one light quark and evaluated near the 
mass-shell of the heavy quark, turned out to be by almost an order of magnitude larger 
than for light quarks Q]. Thus the actual size of radiative corrections in nonleptonic 
decays is essentially not fixed and deserves a detailed study. We address this problem in 
the present paper. 

Apart from the pursuit of completeness of formulas, our study has been fuelled by 
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an intriguing conflict of tlieory with experiment, concerning tlie semileptonic brancliing 
ratio of the B meson: over the last few years the measured semileptonic branching ratio 
T0| has consistently turned out to be (15-20)% smaller than theoretical expectations, and 



the reason for this discrepancy is not understood. This problem has attracted already 
considerable attention among theorists |TT|, |T2|, |T3|, [1^. We believe that the knowledge 



of radiative corrections to nonleptonic decays with full account for the c quark mass is 
imperative for a quantitative understanding of the problem. 

A principal new result of this paper is the calculation of the QCD radiative correction 
to the inclusive decay mode 

• 6 — i> cud 

including full dependence on the c quark mass and the renormalization group improvement 
to the two-loop accuracy. This decay dominates the hadronic B meson width and is crucial 
for understanding the semileptonic branching ratio. We complete the calculations of by 
studying relevant corrections due to operator-mixing and obtain analytic expressions for 
all radiative corrections with one massive quark in the final state, which so far have partly 
been only available in form of one-dimensional parametric integrals |^. As a by-product 
of this analysis, we obtain an analytic expression for radiative corrections to the decay 

• b ^ UTV 

taking into account the r lepton mass. Decays involving the r leptons are potentially very 
interesting, since they may be used to constrain certain extensions of the Standard Model 



and attract an ever increasing experimental and theoretical interest ||T5 
In addition, we estimate QCD corrections to the decay 



considering the colour-favoured contribution, and argue that in this case the radiative 
correction is strongly enhanced when the c quark mass is included. When appropriate, we 
compare our results with those of Hokim and Pham. 

With realistic values of the b and c quark masses, the nonleptonic rate h — > cud is 
increased by approximately (4-8)%, depending on the renormalization scale, if a non-zero 
c quark mass is taken into account which has to be compared with a (4-7)% increase for the 
semileptonic h — >■ cev rate. Thus, the ratio T{h ^ cud)/T{h ecu) is only weakly affected 
by taking into account the c quark mass in the radiative corrections. This compensation 
is due to the specific structure of the effective weak Lagrangian, which supresses the gluon 
exchange between the c and the light quarks. The r lepton mass corrections to the decay 
b —>■ UTU turn out to be of the same order of magnitude as c quark mass corrections to the 
semileptonic b — > ceu decay. 

Although the calculations are straightforward in principle, in practice they turn out to 
be rather tedious and involve several delicate points which we discuss in detail below. We 
use a conservative technique and calculate total imaginary parts by means of Cutkosky 
rules, i.e. we sum up contributions corresponding to real and virtual gluon emission, and 
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introduce small masses for the gluon and the light quark to avoid infra-red (IR) divergences 
due to the presence of massless particles in the final state. The virtual corrections are 
calculated using dimensional regularization, a framework that requires special care with 
the definition of the 75 matrix. We use the so-called naive dimensional regularization 
(NDR) scheme, in which 75 anticommutes in D dimensions. This scheme has been shown 
to yield correct results for the two-loop anomalous dimensions in the relevant terms of the 
effective weak Lagrangian . In the limit of zero masses of all final state particles our 
results coincide with those of Altarelli et al., |0, obtained in the dimensional reduction 
scheme, and with those of Buchalla, obtained in the 't Hooft-Veltman scheme. 

Our paper is organized as follows: in Sec. 2 we recall the general framework for the calcu- 
lation of weak processes to next-to-leading order in QCD, list the necessary contributions 
to be calculated, explain the particular calculation procedure which we use throughout this 
paper, and discuss various subtle points. Sec. 3 contains our main results: the analytic 
expressions for the radiative corrections for non-zero c quark (r lepton) mass, along with 
simple parametrizations. In Sec. 4 we estimate the effect of finite mass corrections on the 
decay h ccs. Sec. 5 contains the numerical analysis and a discussion of various phe- 
nomenological applications. Sec. 6 completes the paper with a summary and conclusions. 
Technical details of the calculation are given in the appendices. 



2 Theoretical Framework 

This section is mainly introductory. We collect the necessary definitions and describe the 
main steps in the calculation of inclusive decay widths to next-to-leading order. 

2.1 The Effective Lagrangian 

Let us consider a Standard Model nonleptonic weak process induced by charged currents. 
In the limit of infinite W boson mass mw and to leading order in the weak coupling, the 
charm-changing nonleptonic decays of beautiful particles can be described by the effective 
Lagrangian 

- ^ {CiiMif^) + c2{Mm ■ (2.1) 

Here O2 is the colour-singlet operator 

O2 = {du)v^A {cb)v-A (2.2) 

appearing in the tree-level Lagrangian, whereas the non-singlet operator {a and f3 are 
colour-indices) 

Oi = {daUi3)v-A {Ci3ba)v-A (2.3) 
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is generated by the strong interaction. The operators (9^ are obtained from the Oi by the 
replacement d ^ s. The Wilson-coefficients Ci{^) can be calculated at the scale /i = my/ 
and then evolved down to scales <C mw by solving the corresponding renormalization 
group equations. Since the operators Oi and O2 mix under renormalization, it is convenient 
to introduce the operators 0± = ((!?2±C^i)/2 which renormalize multiplicatively.Q To two- 
loop accuracy, the corresponding Wilson-coefficients are given by [^, |16| 

C±(/i) =C2(/i)±Ci(/i) 



Here 



^±(/i) 



, .-|7f/(2/3o) 



(2.5) 



is the solution of the renormalization group equation for C± to leading logarithmic order, 
whereas the term proportional to as{mw) — C(s{n) takes into account the renormalization 
group improvement to next-to-leading order. The 7^!'' are the coefficients of the anomalous 
dimensions of the operators 0± and 0'_^, 

,(0)^ , a) {^Y 



with [|T6 



7± = 7r^+7i"(^^j (2.6) 

7f = 4, 7!" = -8, 7i" = -7 + i n„ ^i" = -14 - ? (2,7) 

in naive dimensional regularization with anticommuting 75 and nj running flavours (five 
in our case). The jSi are the ffist two coefficients of the QCD /3-function, 



/?o = ll-^%, /3i = 102-yn;. (2.8) 

Finally, the B± are the matching coefficients ensuring that up to terms of order al{mw), 
matrix elements of the effective Lagrangian calculated at the scale /z = myr equal the 
corresponding matrix elements calculated with the full standard model Lagrangian. Fierz- 
symmetry implies 

S.^±S^. (2.9) 



^Hereafter we imply that the renormahzation scheme obeys Fierz-symmetry. Within dimensional reg- 
ularization this requires certain restrictions on the particular choice of the so-called evanescent operators, 
see raiTTllTsl. 
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where iVc = 3 is the number of colours. In the NDR scheme one finds 



|I6| 



5 = 11. (2.10) 

Note that both the matching coefficients and the two-loop anomalous dimensions depend 
on the renormalization scheme. It proves, however, feasible to combine them into the 
scheme-independent quantity 0, |^, [1^ 



R 




■± 



10863 - 1278n/ + SOnf 15021 - 1530n/ + 80n^ 

Ra- — ; 77! 1 R— — ; tt: 1 (2.11) 

+ 6(33-2n/)2 ' 3(33-2n/)2 ' ^ ' 

in terms of which the Wilson-coefficients read 

C^{^^) = L±(/.) |l + R^ + ^^^bA. (2.12) 

In this form, all the remaining scheme-dependence is absorbed in B± and is to be cancelled 
by the scheme-dependence of matrix elements of the corresponding operators. 



2.2 The Inclusive Decay Rate 



The total inclusive decay rate of a B meson into hadrons Xc, corresponding to the ele- 
mentary subprocess b —>■ cud + cus, can be expressed as the imaginary part of a transition 
operator: 



r(s ^ X, 



TUB 



rriB 



En 



i=l 



(27r)32Ei 



lmiJd'^x{B\T[C 



AC= 
cfT 



lit, 



X 



\0)\B) 



(2.13) 



where the sum runs over all possible final states. Starting from this expression, it can be 
shown that in the limit of infinite b quark mass the B meson decay rate coincides with 
the decay rate of a free b quark, which can be calculated perturbatively. Recently, much 
attention has been paid to understand the structure of nonperturbative corrections to that 
simple picture which are suppressed by powers of the b quark mass . Yet in this paper 
we mainly address the question of how to calculate radiative corrections to the quark level 
decay and thus we proceed in a purely perturbative framework. 

Summing over the squares of the relevant elements of the CKM matrix and making use 
of its approximate unitarity in the first two generations, the hadronic decay width is equal 
to the width of the semileptonic decay b —>■ ceu, 



mr 



1927r3 



iK^lViKK), 



(2.14) 
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multiplied by a colour factor 3 and up to radiative corrections. Here we have neglected the 
light quark masses and the strange quark mass and introduced the notation /i(a) for the 
tree-level phase-space factor: 

/i(a) = 1 - 8a + 8a^ - - 120^ In a. (2.15) 

By the quark masses mb and rric we hereafter understand the pole masses. 

The expression for the inclusive decay width becomes more complicated if one turns on 
the renormalization group running of the Wilson-coefficients in the effective Lagrangian 
and calculates radiative corrections. To next-to-leading order accuracy one has to take into 
account the set of diagrams shown in Fig. |l|. For practical calculations it proves convenient 
to split off the colour algebra and trivial multiplicative factors. Then the operators Oi and 
O2 become identical and the four-fermion vertices in Fig. |l] can be interpreted as insertions 
of O = {du)v-A{cb)v-A or {su)v-A{cb)v-A where no summation over colour is implied. 
For definiteness, the diagrams in Fig. |l] represent the corresponding contribution of 

spin-averaged 

(2.16) 

to the forward-scattering amplitude for an on-shell b quark. The colour-factors of each 
diagram are given in Tab. |l]. 

Taking everything together, we can express the decay rate for b cud + cus as 



r(fe — >• cud + cus) = Fo 



+ ^{2L^(/.)5+ + L2_(^)i?_} 

+ -ALM + L.m'l^{Ga + G,} 
4 o vr 

4 3 TT 



(2.17) 



where Ga, G^, Gc, Gd and G^ are functions of the ratio rric/mh, and are obtained by adding 
the contributions of particular diagrams, which we label by Roman numbers: 



TT 

a. 



Ga = Kim (II + IF + III + IX + IX^ + XII), 



Gb = Klm{W + V + VII), 

TT 
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a 



Gc = Kim (II + II'^ + V + XI + XI^ + XII), 

TT 



as 



— Gd = Kim (III + IV + VI), 



TT 



— Ge = Kim (VI + VIII + X + X^ + XI + Xl^). (2.18) 

TT 

Here 

1 Q27r^ 

K= 2/ 2v (2-19) 



Note that the colour-factors are split off and shown explicitly in (|2.17| ) 



It is easy to verify that Ga and Gb are manifestly gauge-invariant and do not depend on 
the scheme used for 75 . The same is true for Gc and Gd provided the renormalization scheme 
preserves Fierz-symmetry (cf. the next subsection). Taking into account the relation (|2.9| ), 
the contribution of matching coefficients in the second line in (|2.17|) can equivalently be 
rewritten replacing Ge in the last line hj Ge + B, which is again a scheme-independent 
quantity. 

The reason for splitting the full radiative correction to the nonleptonic inclusive width 
into several pieces is that some of these expressions also appear in the semileptonic decay 
widths. For example, the b — >■ ceu decay rate can be written in terms of the function Ga 
as 

Tib ^ cel.) = To |l + ^ ^ G„| (2.20) 

and, owing to Fierz-identities, the function Gc determines the radiative corrections to the 
decay b utu, 

T{b - uTu) = To |l + ^ ^ a} , (2.21) 

with the obvious replacements Vcb — > V^b and a = {mc/mbY {mT-/mbY in To and G^- 
Since the radiative correction to the semileptonic decay b — >■ cev is known analytically for 
arbitrary c quark mass, we can incorporate the answer given in [Q directly in the expression 
for the nonleptonic width and do not need to calculate the function Ga anew. 

Finally, the calculation of Gb and Gd is simplified by using the known results for the 
radiative correction to the vector correlation function with one massive quark: 

n^, = ijd^x e*^" ( I r[g7^Q(x) Q7.g(0)] | ) 

poo , 

= / o {'l^^Ql^Pl (s) + QiJivPXis) \ + constants . (2.22) 

Jml S — (1 

The expressions for the spectral densities (s) can be obtained from we have also 
calculated them directly using Cutkosky rules (cf. App. C). Thus, essentially only Gc and 
Gc remain to be calculated, i.e. the diagrams II, V, VI, VIII, X, XI and XII. 
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It is worthwhile to note that the expressions for Ga and Gb are available in form of one- 
dimensional parametric integrals for arbitrary masses of all the three final state quarks from 
00, corresponding to the upper and lower vertex corrections, respectively, in the notation 
of this paper. Answers for Gc and Gd can be inferred from the expressions given in [0] by 
Fierz-transformations. However, the calculation of Ge involves different Feynman diagrams 
than those calculated in and is new. 



2.3 Calculational Tools 

In this paper we use the well-known dispersion technique and calculate the imaginary 
parts of the necessary diagrams by explicitly summing the contributions of different dis- 
continuities, corresponding to real gluon emission and virtual corrections. Infra-red (IR) 
divergences, which are bound to arise due to the presence of massless particles in the final 
state, are regulated by introducing non-zero gluon and light quark masses. According to 
the famous Bloch-Nordsieck theorem, all IR divergences cancel in the sum of all possible 
discontinuities for a given diagram, so that in the final answer for the total imaginary 
part one can put the regulator masses to zero. In principle, one could escape the problem 
with the IR divergences altogether calculating the matrix element in question in Euclidian 
space, continuing analytically to Minkowki space and taking the imaginary part. This 
technique is very powerful for massless quarks, but, as we have found, is not very useful 
when quark masses have to be taken into account. As for diagram V, we prefer to express 
the subdiagram containing the self-energy by a regularized dispersion relation and only 
then compute the discontinuity. 

We further use dimensional regularization and the MS subtraction scheme to deal 
with ultraviolet divergences, which give rise to the renormalization of operators and quark 
masses. 

Due to the presence of 75, the use of dimensional regularization involves subtleties that 
have received a lot of attention over the last years. We have found that for our purposes 
it is sufficient to use the simplest prescription with anticommuting 75 in D dimensions, 
NDR. It is known that this scheme can yield algebraic inconsistencies when the amplitude 
in question contains closed fermion loops. Our method of calculation explicitly avoids 
computing such diagrams by the appropriate use of Fierz-symmetry. Indeed, by applying 
Fierz-transformations one can always achieve that any diagram contains no dangerous 
closed fermion loops,^ as illustrated in Fig. |^. Subtleties are bound to arise, however, if 
the Fierz-symmetry is broken by regularization. Thus, we use the results of Buras and 
Weisz [0 who have checked that in the NDR scheme the renormalized operators Oi and 



^We use the expressions given in the second of Rcfs. ||^ and imply that the factors In A/w appearing 
in Eq. (3.35) of this paper should be understood as (In A)/w, as we have found by comparison to the first 
of Rcfs. 0, and to the numerical results given in 

■^In fact, closed fermion loops with only one external momentum and at most two open Lorentz-indices 
do not pose a problem, since Tr(7^ 7^/75) vanishes in any dimension. 
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O2 have the same properties under Fierz-transformations as the bare ones, 



02^0^{c^d), Ci^C2(c^rf). (2.23) 



Furthermore, with the particular choice of evanescent operators made in |T6| this property 
is vahd diagram by diagram: after subtraction of counterterms, renormahzed Feynman 
diagrams that can be related to each other by means of Fierz-transformations in four 
dimensions coincide identically. Note that we are only interested in the imaginary parts 
of the multi-loop diagrams of Fig. |I|, which are given by the product of at most one- 
loop and tree-level amplitudes, followed by phase-space integration. Thus, evaluating 
the product with renormalized amplitudes, one can avoid the contributions of evanescent 
operators altogether, provided their choice is consistent with the calculation of the one- 



loop matching coefficients B± by a naive projection onto four dimensions, see for 
details. We have checked that the alternative procedure for calculating the whole three- 
loop diagram in dimension D, taking into account the counterterms due to the evanescent 
operators as given in [jl6|, yields the same result. Last but not least, we have found that 
in the limit rric 0, our answers coincide with those obtained by Altarelli et al., [Q, 
using dimensional reduction, and by Buchalla, 0, using the 't Hooft-Veltman scheme, cf. 
App. |B.8| . More details about the calculation and some useful intermediate formulas can 
be found in appendices A to C. 



3 Analytic results 

In this section we give the explicit expressions for the next-to-leading order corrections 
to the nonleptonic rate b — > cud, the functions Gx defined in Eq. ( p.l8| ). The answers 



for all necessary single diagrams are given in App. B. Throughout this section, we use 
a = {nic/mbY; the factor K was defined in (|2.19| ) and the phase-space function /i in 
( p.l5| ). All results are given in terms of pole masses. 

As mentioned before, Ga is proportional to the semileptonic width r{b ^ cev). From 
the analytical result given in P] we obtain: 

l?-^'^ = -St^^ {-(1 - - 956a + 75a^) + 4(1 - a') 

X (17 - 64a + 17a^) ln(l - a) + 4a(60 + 270a - 4a^ + I7a^) Ina 

- 48(1 + SOa^ + a^) In a ln(l - a) + 12a2(36 + a^) In^ a 

- 384a3/2(l + a) [vr^ - 2 In Ina + 4L2 f-a^/^A _ / i/2\ 

[ 1 + a^l^ ^ ^ ^ 

- 12(1 + IGa^ + a^)(-7r2 + 6L2(a))| . (3.1) 
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We next calculate Gh and Gd by exploiting the result for the vector correlation func- 



tion given in ||19[. To get G^, we first apply a Fierz-transformation to both four-quark 
vertices of the diagrams III, IV and VI, cf. Fig. ^ We then express the closed fermion 
loops by twice the radiative corrections to the vector correlation function^ II^i. defined 
in Eq. ( p.22| ) in their dispersion representation. This effectively replaces the loop by the 
propagator of a pseudoparticle, whose mass is just the invariant mass of the original two 
fermions, multiplied by the weight functions p\ given in App. C and followed by integra- 
tion over phase-space. By cutting the two remaining propagator lines we obtain the sum 
of imaginary parts of the diagrams III -|- IV + VI: 



K TT 



G, 



ds{mi-sf mlpf\s)~2p'^\s) 



(3.2) 



The integral can easily be done with the formulas given in App. A. The final result is: 

1 a. 



1? — 

K TT 



2 fi 

92167r5 



{(l-a)(l 



8 - 476a - 539a^ + 75a'' 



+ 4a2(36 + 8a-a 



2^ 
1 I 



3 In^ a 



- 2(1 - c?) [31 - 320a + 31a^ - 12(1 - 8a + a^) In a] ln(l - a) 

- 2a(132 + 90a - 308a^ + 31a^) In a 

-f 24(2 - 16a - 36a2 + 8a^ - a^ 12a2lna)L2(a) 
+ 864a2[C(3)-L3(a)]}. 



(3.3) 



For Ghi the fermions in the loop are massless, so that 



K TT 



G. 



Anm^ Jo 
— 2{mi, + — s)p2 [s 



2 2\2 / 2 , 2\ 



massless 



massless 



(3.4) 



2 fi 

5127r5 
where w is defined as 



/i(a) 



w(x, y, z) = (x^ + y'^ + z'^ — 2xy — 2xz — 2yz) 



1/2 



(3.5) 



Since K is proportional to l//i(a), the above result simply means that Gb is a constant 
independent of a. 

"'Note that the vector correlation function equals the axial vector correlation function in perturbative 
QCD. 
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Summing up the diagrams II, 11^, V, XI, XI^ and XII one gets 
1 9s^b 



kV'^' = frn^S - - 539. - 476a^ + 18.= 



- 4(3 - 24a - 36a^ + 16a^ - 2a^)7i^ 

- 2(1 - a^) [31 - 320a + 31a2 + 12(1 - 8a + a^) In a] ln(l - a) 

- 2a(12 + 90a - 188a2 + 31a^ - 48a vr^) In a 

- 24(1 - 8a + 36a2 + 16a^ - 2a^ + 12a2 In a)L2(a) 

+ 864a2[L3(a) -C(3)]}. (3.6) 

Finally, summing up the diagrams VI, VIII, X, X^, XI, XI^, we obtain the scheme- 
independent quantity + B: 

^— {Ge + B) = S^S^ - «')(41 - 488a + 41a2) - 16a(l + 3a + a^)^^ - 6L2(a)] 
K 71 15367r^ 

- 8(1 - a^)(l + 28a + a^) ln(l - a) + 4a2(120 - 32a - Sa^) Ina 

+ 48a^ln2a + 2(61n(m2//i2) + ll)/i(a)}. (3.7) 

Note that Ge contains a logarithm of /x, namely Gg = Q ln(m^//i^) + . . ., which just cancels 
the remaining fi dependence of the leading order term in Eq. ( |2.17|) , as shown by an explicit 
expansion around = m?,: 



2L^(/x) + L2_(/x) = 2Ll{m,) + L'_{m,) - ^ {L^(m,) - Ll{m,)} | 6 In ^. (3.8) 



This completes the list of entries required to evaluate the radiative corrections to the 
nonleptonic decay b — >• cud. We have checked that our expressions for Gc and Gd agree 
numerically to the results given in As mentioned above, the answer for Ge cannot be 
inferred from this paper. 

For practical calculations, it may be useful to give approximations to the above formulas 
for small a: 

31 

{Ga + Gb)fi{a) = — - TT^ - a[80 + 24 In a] + 32a^/2 tt^ - a2[273 + IGtt^ - 18 In a 

+ 36 In^ a] + 32a^/2 ~ g [^^^ - 57 Ina] + O (a^/^) , (3.9) 
31 

{Ge + Gd)fi{a) = — -7r2-8a[10-7r2 + 31na] - a2[117 - 2471^ 
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+ (30 - Stt^) In a + 36 In^ a]--a^ [79 + 2n^ 



621na + 61n2a] + C(a^) 
51 



mt 



{Ga + Gt + Ge + B)Ma) = 61n^ + ll /i(a 



7r^ + 8a[21-7r^-31na] 



+ 32a^/2 - a^iUl + AOtt^ - 258 In a + 36 W a] 



(3.10) 



+ 32a^/2^ 



^ -a^[305 + 187r2 + 301na-541n2a] 
9 



+ 0[a 



7/2 



(3.11) 



These approximations are valid for a < 0.15, which covers a reahstic range of values 

The corresponding approximate formula for the semileptonic decay b cev is given 
in p[. We repeat it here for completeness, converting to our notations, cf. ( p.20| ): 



25 

Gah{a) = — - TT^ - a[68 + 24 In a] + 32a^/V^ - a^llQ-K"^ + 273 - 36 In a + 36 In^ ( 



+ 32a^/V2 



1052 152 



In a 



9 3 

Similarly, we obtain for h uru, cf. ( |2.21| ): 
25 

GJi{a) = —-71^ - 8a[Q - tt^] - 20^(15 - Qtt^ + (15 - An^) Ina + 36C(3)] 



(3.12) 



+ a-" 



20 



167r2 



(3.13) 



where here a has to be interpreted as {rrir/mby. 

These expressions complete our formulas of radiative corrections to b quark decays with 
one massive particle in the final state. 



4 Estimate of Finite Mass Corrections to r(6 



CCS 



The structure of the decay rate r(6 — >• ccs) is very similar to r(6 — >■ cud), except for the 
additional contributions of penguin diagrams. Since penguin operators enter the effective 
Lagrangian only by mixing, their Wilson-coefficients are small and thus penguin contri- 



butions can usually be neglected in tree-level decays, cf. |11|. Defining functions in a 
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completely analogous way to the Gx that enter V{h cud), Eq. ( p.l7| ), the decay rate of 
b CCS + ccd can be written as 



CCS + ccd) = T'q 



2L^(/i) + L^(/i) + 



as{mw) - ctsifJ') 



271 



{2Ll{f,)R^ + L'_{f,)R_} 



+ 



27r 



4 O 71 

4 6 71 



+ penguins. 



(4.1) 



The tree-level decay rate Fg is defined as in Eq. ( p.l4| ) except for the phase-space factor. 
Neglecting the strange quark mass, /i(a) has to be replaced by 



/i(a,a) = y/T 



T^{l-Ua-2a'-12a')+24a^ (l - a^) In ( |^^ |_^^ ) 



(4.2) 



with a = {rric/mb)^. With = 0.2 GeV, (|4.2| ) is accurate to better than 5% for a < 0.35^. 
In the subsequent estimate of r(6 ^ ccs), however, we take into account a finite rris and 
use the three-particle phase-space factors calculated in [pO[] . 

Ha parametrizes the corrections to the decay rate b —>■ ctv with the "r" and the c quark 
degenerate in mass, whereas He now describes the corrections to a decay of type b utu, 
where the light quark u is still massless, but the leptons have equal finite masses. 

The calculation of Hj, requires only expressions already obtained in this paper. From 
Eq. (|3.4|), where Gb is expressed as a single integral over spectral functions, we immediately 
obtain: 



K' 71 



4:7iml 



(mi,-mc)^ 
-,2 



ds { [{ml - mlf - {ml + m^)^] ^ {s) 

- 2{ml + ml- s)pf{s)]w{s,mlml) (4.3) 

with w defined in ( p.5|) and pf^ given in App. C. K' equals K, Eq. ( |2.19|) , except for the 
obvious replacement of phase-space factors, /i(a) — > f\{a,a). Hj, is extremely sensitive to 
the value of the c quark mass as clearly visible from Tab. § recall that the corresponding 
quantity for b — > cud, Gb, was constant in m-c! The numerical values given in Tab. ^ agree 
with the corresponding ones obtainable from 0. In the table we also give the values of Ha, 
He and H^ as functions of mc/ mb calculated from the formulas in [0 . He cannot be obtained 



13 - 



from the work of Hokim and Pham, so for the present estimate we insert ~ Ge, Eq. 
( p.7|) , and estimate the error introduced by this procedure by twice the difference between 
Ge{a) and Ge(0). From Tab. |^ it is clear that both Ha + Hj, and He + Hd are extremely 
sensitive to rric. For the reahstic value rric/mb ~ 0.3, both Ha + Hf, and Hc + Ha are positive 
and by a factor four to five larger than the corresponding Gx- In Tab. ^ we also give the 
decay rate r(6 ccs) including the error in ifg- For nT-c/^b = 0.3, the finite c quark mass 
effects in the radiative corrections increase the decay rate by 30%. 



5 Numerical Results 

Let us now turn to the numerical analysis of the results obtained in Sees. 3 and 4. To 
do this, we first rewrite the decay width Eq. (|2.17|) in a form that is more convenient in 



the present context (a = (mc/mf,)^). Introducing scheme-independent coefficients Cij by 
cii = Gc + Gd, Ci2 = Ga + Gb + Ge + B and C22 = Ga + Gb, we have 



r(6 — > cud + cus) = Tq 



+ ^ {L4f^) + I ^ C22{a) + ^ - I ^ cn{a) 



(5.1) 



3 vr 

= 3ror/(/i)J(a,/i). (5.2) 
Here we use the same notations as [|11[] with 

vi^,) = ^-{2Lli^,) + Lli^,)}, (5.3) 

whereas J = 1 + 0(q!s) covers the next-to-leading order terms. 

In Fig. ^, the coefficients cn, Ci2(/i = mb) and C22 are plotted as functions of rric/mb] we 
have also tabulated the numerical values in Tab. ^ A range of realistic values of rric/mb, 
0.26 < rric/mb < 0.32, is indicated by a grey bar. The figure shows clearly that the colour- 
favoured coefficient C22 is not very sensitive to the c quark mass. This is due to the fact that 
Gb = 3/2 is constant in rric, whereas the term Ga, which also determines the semileptonic 
decay width, exhibits an only mild dependence on m^. The coefficient C12 increases more 
steeply in rric, but is suppressed in the rate by a small Wilson-coefficient. Both C12 and C22 
are finite for rrib. This is different for Cn: this function diverges like ~ ln(mfe — m^) 

as shown in Tab. ^. Note, however, that the decay rate still vanishes for rric = ^b- For cn 
and C22 we found agreement with the integrals given in |^. 

We next investigate the behaviour of r(6 cud) as compared to r{b ecu). To this 
purpose, we introduce 

/(a,/i) = l + ^^G,(a), (5.4) 

3 TT 
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so that r(6 ceu) = ToI{a,fi). In Fig. Ka) we plot I{a, fi)/I{0, fi) as function of 
mc/mf, for three different values of the renormalization scale /i. Fig. ^(b) shows the ratio 
J(a, /i)/J(0, /i) as function of rric/mb. For rric/mf, ^ 0.3, we find an increase of the nonlep- 
tonic decay rate b — > cud by (4-8)%, depending on the renormalization scale. Comparing 
the two plots, we find that in the range of realistic values of rric/mb, which is emphasized 
by a grey bar, the enhancement by finite mass corrections is approximately equal for semi- 
and nonleptonic decays, whereas for mc/mh near one, the logarithmic divergence in cn 
clearly dominates. 

Whereas up to now we have discussed only quantities that are defined on a purely 
partonic level, one clearly is more interested in predictions of hadronic decay rates and 
branching ratios. Though in general this requires the knowledge of hadronic matrix ele- 
ments, matters simplify considerably for inclusive decays of heavy hadrons. Without going 
into further details, we just mention that the theory of inclusive decays of heavy hadrons 
has experienced considerable progress over the last years when another field of applica- 
tion of the well-known operator product expansion (OPE) technique as an expansion in 
inverse powers of the heavy quark mass was opened, see [|T], ^ . The situation is even more 
favourable than in other fields since the leading term in the heavy quark expansion does 
not involve any hadronic uncertainties, but simply coincides with the underlying free quark 
decay process. Hadronic corrections enter only at second and higher order in the expansion 
with a natural size of ~ 1 GeV^/m^ ~ 5%. 

To second order in the heavy quark expansion, these corrections can be expressed in 
terms of two matrix elements, 

2mB\i = {B\K{iD)%^\B), 

em^As = ( 5 I a^^F^'K \B), (5.5) 

where is defined as = e^^*''"^h{x), b{x) being the b quark field in full QCD, is the 
four- velocity of the B meson, its mass and F'^" the gluonic field-strength tensor. 
Whereas A2 is directly related to the observable spectrum of beautiful mesons, 

A2 ~ ^ (m|. -m|) = 0.12GeV^ (5.6) 

the quantity Ai is difficult to measure, cf. |^. Physically, — Ai/(2mf,) is just the average 
kinetic energy of the b quark inside the meson. At present, only a QCD sum rule estimate 
is available, according to which Ai ~ —0.6 GeV^ [22|. This result has been met with 



caution (see, e.g. ||2^), since it corresponds in fact to a surprisingly large momentum of the 
b quark inside the meson of order (700-800) MeV. However, in a recent series of papers, cf. 
pl| , p^ , an upper bound on Ai has been derived, to wit Ai < —0.4 GeV^, which appears to 
be in nice agreement with the QCD sum rule prediction. In our analysis we thus conform 
to the value Ai = -(0.6 ± 0.1) GeV^. 

Taking into account hadronic corrections, the nonleptonic decay rate of a B meson into 
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a single charmed hadronic state can be written as 



STo r]{fi)J{a,fi) \ 1 + 



Ai 3A2 6(1 - a)^ A2 



2ml 



+ 



2ml 



mf 



4(1 -a)3 A2 
/i(a) ml 



a. 



\ mt mi 



(5.7) 



For the semileptonic rate one obtains: 



T{B ^ X,eu) = To 



/(a,^){l + A^ + 



3A2 
2ml 



6(1 -g)^ A2 



mt 



\ s 1 2 ' S 



(5.^ 



In evaluating these decay rates, it is crucial to minimize the number of independent 
parameters. In doing so, we take advantage of the fact that the difference between heavy 
quark masses is fixed in the framework of the heavy quark expansion: 



mb — mc = ms — mo + 



Ai + 3A2 











--) 






mc/ 





The only quantity remaining to be fixed is then m^ or mc 
spectroscopy and choose the most conservative range| 



With these values we find: 



4.5 GeV <mb< 5.1 GeV. 

T(B ^ Xc) 
^ ^' 4.0 ±0.4. 



(5.9) 

We prefer to take m^ from 

(5.10) 

(5.11) 



This result is nearly independent of the c quark mass; the error is entirely due to the 
dependence on the renormalization scale, mb/2 < fi < 2mb- Note that it is mainly the 
mass-independent Gb that enters the radiative corrections, whereas Ga cancels, which 
explains the small sensitivity of ( ^.11 ) to mc- 

Finally we calculate the semileptonic branching ratio as 



B{B Xev) 



T{B Xeu) 



E^=e,M,r r(5 ^ Xeue) + T{B -.Xc) + T{B^ Xc 



(5.12) 



Previous analyses yielded B{B Xeu) > 12.5% |TT], |T^ which is considerably bigger than 
what is measured experimentally, B{B Xeu) = (10.43 ± 0.24)% [^. Yet these studies 
were lacking the finite quark mass effects on J {a, n), the radiative corrections to T{b ccs) 
and the corrections to r[b — > crz/). Taking them into account, we find B{B Xeu) = 
(11.6 ± 1.8)% where the error comes from the uncertainty in the renormalization scale, in 



Since we give a broad range of values, we are not concerned with the intrinsic uncertainty in the 
definition of the pole mass which is caused by renormalons and estimated to be of order (50-200) MeV, 

p5i. 
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Ai, in the quark mass and in the value of as{mz), where we used as{mz) = (0.117±0.007) 
The functional dependence of the semileptonic branching ratio on rrih with fixed by 
j9[) is also plotted in Fig. |^. Note that in r(_B ^ Xcc) in Eq. (|5.12| ) we use the full tree- 
level phase-space factor calculated in pOf including a strange quark mass rris = 0.2 GeV. 
With that value of m^, the rate is decreased by typically 5%, which is a much smaller effect 
than taking into account the c quark mass in the radiative corrections. 

In Tab. ^ we compare the semileptonic branching ratios obtained in |jTl], |l3l with ours, 
using the same input parameters. It turns out that, although the introduction of the heavy 
quark expansion by Bigi et al. leads to a reduction of the branching ratio by 0.3% compared 
to the quark level analysis of Altarelli and Petrarca, the finite c quark mass effects in the 
radiative corrections to the nonleptonic widths are by more than a factor three bigger 
and yield an additional 1.0% reduction. This fact shows the paramount importance of 
perturbative corrections in the heavy quark expansion. 

Yet, we consider this analysis as a rather preliminary one. To make more definite 
predictions, one has to reduce the scheme-dependence. We plan to come back to this and 
related questions in a separate publication 



6 Summary and Conclusions 

In this paper we have calculated the radiative corrections to the free quark decay b 
cud as the imaginary part of the relevant forward-scattering amplitude including the full 
dependence on the c quark mass. We have performed the calculation in naive dimensional 
regularization (NDR) with anticommuting 75. The key-observation that allowed us to 
employ this scheme despite of its known deficiencies was that diagrams that are ambiguous 
in NDR can be related to well-defined ones by means of Fierz-transformations. In the 
limit rric — >■ 0, our results agree with known expressions obtained in other schemes. As far 
as comparison is possible, our results also agree with those of 0. 

For realistic values of the b and the c quark mass, we find a moderate increase of about 
(4-8)% of the decay rate r(6 — cud) with respect to the limit of vanishing c quark mass, 
depending on the renormalization scale. 

We also have estimated the increase of r(6 — ^ ccs) by finite mass effects to be about 
30%. The calculation of the full radiative corrections is under study. 

In the framework of the heavy quark expansion of inclusive decays of heavy hadrons, the 
knowledge of finite quark mass corrections to nonleptonic decays is crucial for an improve- 
ment of the theoretical prediction of various measurable quantities like r{B X^) /^{B —>■ 
XcCu) and B{B — > Xeu). We have shown that the ratio of non- to semileptonic decays 
of B mesons to single charmed final states remains nearly unaffected by finite mass ef- 
fects, whereas the semileptonic branching ratio is reduced by 1.0% compared with a recent 



investigation by Bigi et al., |T^. The finite mass terms in the radiative corrections are 
thus of the same importance as the hadronic corrections introduced in [Q, Q and severely 
reduce the scope of possible new physics in nonleptonic B decays. Yet, the analysis of the 
semileptonic branching ratio involves delicate points like the choice of quark masses and 
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estimates of higher order radiative corrections; we will address these points in a separate 
paper p8 |. 
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Appendix A 

In this appendix we collect some intermediate formulas which can be useful in reproducing 
our results and related calculations. 



A.l One— Loop Integrals 

We define the relevant one-loop integrals as 

(K K K ) = f (1, kp, kpK) 

yjx,jxp,jxp,) J ^27r)/3[(A;-pi)2_M2][(A; + p2)='-/i='][A;2-A2]' ^ ^ ^ 

where D is the space-time dimension. Throughout the appendixes, it is always assumed 
that pI = ^ pI = fi^ ':$> . M (/i) plays the role of the b or c (u or d) quark mass, 
whereas A is the gluon mass which along with fi is used to regularize possible IR singularities 
of the integrals. The renormalization scale appearing in regularized UV divergent integrals 
is denoted by u. 

The one-loop integrals can be parametrized as 

K=-j^K{M\q'), q = Pi+P2, (A.2) 
Kp = -j^ {-25(M2, q')qp - 2B{M\ q^)p,p] , (A.3) 

Kpa = [A{M\ q')qpq„ + B{M\ q^) [p^pq, + P2M + B{M\ q')p2pP2. 

+ C{M\q')gp^} (A.4) 
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with 

~ 1 f / 



C(M^ = -il^ I ^ + In 4. - 7. - In ^ 



(A.8) 



Only C{M^,q^) is UV divergent. Both K{M^,q^) and B{]VP,q^) are IR divergent, while 
A(M2, ^2) and S(M2, are UV and IR finite. 

For g2 > 71^2 ^j^ggg functions become complex. It then proves convenient to separate 
real and imaginary parts explicitly: 

~ 1 f / ;/2 A/f2 \ 

q^ , 



' > M^) - ^ {l + (l - f ) (in ^ + ) , (A.11) 
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1 / \ 
i3(A/^ > M-') = [in + i^j , (A.12) 

B(M^ ,^ > m = -^p^ {in ^ + 1„ ^ + 2..) . (A.13) 

C■(M^ > M=) = I + 1,1 4^ - 7E - 1.1 

A. 2 Phase— Space Integrals 

Here we collect only the most involved three-particle phase-space integrals. The simpler 
ones, as well as formulas for n-particle phase-space integrals can be found in . 
We define 



/lIPS(p„p„A:) = /|||£||5V-J>i-P2-A:) (A.15) 
with pI = fi^, k"^ = and keep only logarithms of and /x^. 



^ {2(P= -p;)=fteB(j';. P') - ^(rf, P=)} . (A.19) 



A2 + 2pi ■ /t 4P2 I P2 

/urora^- -^{(P^-ri)Re/7(rf.P') + ^'(ri.P=)}. (A.21) 
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{X^ + 2pi-k){X^ + 2p2-k) 



7r2 



-pl)^ReB{pl,P^) + P^ -pI}. (A.22) 



2P2 

In addition to the functions defined in the last subsection, we also have used 



X — y X — y 3 \y — x 



In I In - ![! - 2 L2 f , (A.23) 



X X — y 3 \x 
E'{x,y) = Eiy,x) + ^, (A.24) 

X 

H{x,y) = X — y + y\n —. (A. 25) 

y 



Note that ([AlBD , (lAlOl) , (|A:21|) and (^[H) are IR divergent. 



A. 3 Integrals Related to Polylogarithms 

Our approach yields the relevant diagrams in terms of an integral over the invariant mass 
squared of two quarks, either both massless or one massive and the other one massless. 
Using the formulas given in this appendix, the remaining integral can also be done. There 
appear at most polylogarithms, which are defined as 

Ln(y) = / — L„_i(x), Lo = -ln(l-y), (A.26) 

JO X 

and have the following simple series expansion in the unit disc: 

00 ^,fc 



Ln(?/) = E!^' \y\<^- (A-27) 
fc=i ^ 



We need the generic integrals 



/„ = C dx ln(l - x) In x, (A.28) 

J a 

Jn = ^ dxx''L2ix), (A.29) 

Ja 

Kn = C dx x" ln(x - a) In x, (A.30) 

J a 

Ln = f^dxx''u{^ (A.31) 



for integer n. They can be computed with the help of the following recursion formulas : 

n ^ a" In a , , , , , 

Inin > = —— In-i + — {a + 1 - a) In 1 - a)} 

n + 1 n + 1 



+ —^J dx|x" + x"-^(l-a;)ln(l-x)+na;''lnx}, (A.32) 



Jo = 2 - 2a - — - (1 - a)(l -lna)ln(l - a) + a In a + L2(a), (A.33) 
6 

J_i = L2(a)lna + C(3) -L3(a), (A.34) 
In{n< -1) = ^^^|-a"+Mn(l-a)lna-L2(l-a)-iln2a| 

\r|^l°(l-^), (A.35) 

1 Ja X 



n + 



j(n^-l) = ^^1— -a"+^L2(a)+ dx ln(l - a;)| , (A.36) 

n + 1 [ 6 Ja J 

J_i = C(3)-L3(a), (A.37) 



1 da; 



= « / ^ In - + ln(l - x) In a - a"+^/_„_-2, (A.38 



L„ = a"+V_„_2. (A.39) 
Note that the sum in ( |A.35 ) evaluates to zero for n = —2. 



Appendix B 

In this appendix we give some details of the computation of the imaginary parts of the dia- 
grams II, V, VI, VIII, X, XI and XII. We have used Cutkosky rules, hence, we think it may 
be useful for the reader who wishes to reproduce our calculations to give the contributions 
of the different cuts separately. 

Colour-factors are omitted throughout the appendix. In order to get the final answers 
for the diagrams, one has to multiply by the appropriate colour-factors shown in Tab. |l[ 

We have used the NDR scheme with anticommuting 75; for a discussion of this procedure 
cf. Sec. 2.3. We work in Feynman gauge and use MS subtraction. In contrast to App. A, 
the renormalization scale is now denoted by fi. The definition of the b quark mass as pole 



- 22 - 



mass is fixed by the treatment of tlie diagrams II and XII. The c quark mass is defined 
by the treatment of the graphs containing the c quark self-energy; we likewise choose the 
pole mass definition. 

Throughout this appendix, we use a = {rric/mbY. The phase-space factor /i was 
defined in ( p.l5| ). 

B.l Diagram II 

Diagram II is the simplest one. One only needs to multiply the lowest order diagram I by 
the finite part of the b quark's on-shell wave-function renormalization constant, Z2F: 

where T,'{mh) is the derivative of the b quark self-energy with respect to |^ at |^ = nib- The 
UV divergent piece of Z2F can be subtracted from the imaginary part of diagram II before 
performing the phase-space integrals. One finally gets (neglecting the colour-factor Cp): 

The IR divergence in the gluon mass. In X^/ml, is an artifact of the on-shell normalization 
and cancels against a corresponding term in ImXII. 



B.2 Diagram V 



This diagram involves the d quark self-energy. It can be computed in several ways of which 
we sketch only the simplest one. We first observe that the diagrams IV and V are equal, 
since the u and the d quark have equal masses, which immediately allows us to write a 
relation analogous to ( p.2|) for the twice Fierz-transformed diagram: 



ImV 



ds (ml 



[B.3) 



where the pf^'" are the spectral densities of the light quark's self-energy contribution to 
the vector correlation function 11^,^, Eq. ( p.22 ). They are given in Eqs. ( C.8 ) and 
Performing the integral in ([B.3|) then yields: 



ImV = 



2 fi 
184327r' 



; |(1 - a2)(-43 + 296a - 430^) - 480^ 



+ 12/i(a) 



21n(l - a) +ln 



mt 



+ 2880^ In a In(l-a) 



12a^(18 -8a + a^)lna + 2880^12(0)} 



(B.4) 
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B.3 Diagram VI 



Applying Cutkosky rules and using the formulas of App. A, both three-particle cuts yield 
a complex result. In the sum, however, only the real parts contribute, for the cuts are 
complex conjugate to each other. One finds: 



9's 



r"^i ds 



s - mj {rrit, - s) 



{ [2s^ + {ml + ml)s + 2mlml 



X 



[{s-ml) \m.eB{ml,s) -Rei^(m^,s)] + 2ReC(m^, s) + l} 



+ s \2s'^ + {m\ + ml)s — m^ml Re A(ml, s) 



+ s 



2s'^ + {Aml + ml)s + 5mlml Re5(m^,s)|. 



(B.5) 



The pole of the function C{ml,s) in D — A cancels after addition of the counterterm 
diagram. Hence, only the UV finite pieces must be retained. Yet the subtraction of UV 
divergences has to be done with some care: counterterm diagrams must be added before 
performing the phase-space integration, which we explicitly carry out in four dimensions. 
This procedure has the advantage that one can ignore the evanescent operators. Note 
however that the three-particle cuts are IR divergent since they contain the functions 
K{ml. s) and B{ml, s). 

With the help of App. A one can also check that the sum of the two four-particle cuts 



IS 



192n^ml 



ds 



m2 S-" 



— {mb — s) < 2s + {mc + m^)s + 2m^m^ 



X {s-mlf{{s-ml) [Reis:(m^,s) - 2ReS(m^,s)] +E'{ml,s)] 
+ {2ml — ml)s + Am^ml ^(^^1^^) 



— m 



^(s + 2m^)(s + 2m^)ln-^-^ 



- {llmi - 20mi)s^ + m^(19m^ - 22m^)s + SSmjml 



(B.6) 



Obviously, the four-particle cuts are UV finite, because they do not involve any loop 
integration. Note that the IR divergences cancel in the sum of the three- and four-particle 
cuts. The result can be integrated with the help of App. A to yield: 



ImVI 



2 fi 

92 1 67r5 



;i - a) (61 - 165a - 3900^ + 52a^) + 4a^(30 + 8a - 0^)77^ 
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2(1 -a^) 25 - 272a + 25a2-12(l-8a + a2)lna ln(l - a) 



+ 2a(12 + 18a + 2120^ - 19a^) In a - 120^(18 + 8a - a^) In^ a - 12/i(a) In 



+ 



24(2 - 16a - 30a^ + 8a^ - a^ + 12aMna)L2(a) + 864a^[C(3) - Uia)]}. (B.7) 



B.4 Diagram VIII 

Proceeding with diagram VIII along the same hnes, one gets for the sum of the two three- 
particle cuts: 



ImVIII^^-P^''*-) 



f 



ds 



ReA{ml,s) + Re5(m^,s 



+ (s-m^) 2ReB{ml,s) -ReK{ml,s) + 8ReC(m^, s) + - L (B.8) 



while the sum of the two four-particle cuts is 



ImVIII^^-P'^''*) 



9l 



ds 



ml S 



— {ml - sy { {s - mlY 2Re S(m^, s) - Re K{ml, s) 



- {s- mlfE'{ml, s) + {s- ml) — — — - 2mls In * 



C 9 



(B.9) 



Again, UV divergences are cancelled by the counterterm diagrams, whereas the IR diver- 
gences add to zero in the final sum. The latter one can be integrated to give: 



ImVIII 



2 fi 

92167r5 



;i - a) (-220 + 1374a + 1863a^ - 139a^) - 4a^(66 + 32a - a^)?!^ 



+ 2(1 - a) [37 - 379a - 811a^ - lla^ - 12(1 + a)(l - 8a + a^) In a] ln(l - a) 

- 2a(12 - 918a + 596a^ + 35a^) Ina + 12a^(18 + 32a - a^) In^ a 

- 24(2 - 16a - 66a2 - 16a2 - a^ + 12a2 In a) U^a) 



+ 48/i(a) In ^ + 864a2 [U{a) - C(3)] . 
B.5 Diagram X 

With the same notations as above, we get for diagram X the following results: 



(B.IO) 



ImX(3-P--) = /J ^ {ml - s)%s - mlf [2s[A{ml s) + B{ml s)] 



+ 8CK, + i + K - s)[K{ml s)-2B{ml s)] \ . 



(B.ll) 



Note that A{ml,s), B{ml,s), C(m^,s), K{ml,s) and B{ml,s) are now real functions 
because s < . By adding the counterterm diagram one removes the UV divergent piece 
in C{ml, s), as it should be. 
The four-particle cut yields 



+ {ml - sfE{ml s) + {ml - s)H{ml s) - {ml - 5) ^^^""' 



+ s{s + ml) In 



(B.12) 



so that 
Im [X + X^] 



{(1 - a) (-139 + 1863a + 13740^ - 2200^ 



92 167r5 

+ 4(3 - 48a - 660^ + IQcC' - 2a^)n^ - 864a^ [13(0) - C(3)] - 2(1 - a) x 
X fll + 811a + 379a^ - 37a^ - 12(1 + a)(l - 8a + a^) In a] ln(l - a) 



+ 2a(12 + 486a - 236a^ + 13a^ - 48a7r^) In a + 48/i(a) In 

+ 24(1 + 16a + 66a^ + 16a^ - 2a^ + 12a^ In a) L2(a)}. 
Again the result is both UV and IR finite as it should be. 



(B.13) 



B.6 Diagram XI 



ImXI^^-part.) 



™? ds 



ml — sY{s — mlY I 2s^ + {ml + ml)s + 2m^m^ 



X {{ml - s) [2B{ml s) - K{ml s)] - 2C{ml, s) - l} 



2s^ + (m^ + m^i)s — mlml 



A{ml,s) 



— s 



2s'^ + {ml + Aml)s + 5mlml B{ml,s)^ 



(B.14) 
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A{inl, s), B{ml, s), C{ml, s), K{ml, s) and B{ml, s) are real functions and the counterterm 
diagram cancels the pole of C{m1, s) in D — 4. 
The four-particle cut is: 



''i ds 



— (s — m^) < 2s + {m^ + m^)s + 2mjnij {m^ — s) 



X 



{{ml - sf [K{ml s) - 2B{ml s)] + H{ml s) + {ml - s)E{ml, s)} 

12 [22^^ 



+ 2 



^ ml ml — s 



4s^ + {2ml ~ ^&)'^ ~ 2mlml 

+ {8ml - ll^b)^^ - ml{22ml - 7ml)s + 14m^m^] }. 
After integrating over s, the final sum can be written as: 



2 6 

Im [XI + XI^] = liH^ {(1 - a) (52 - 390a - 165a^ + 61a^) 



92 167r5 



4(3 - 24a - 30a^ + 16a^ - 2a^)7r^ 



(B.15) 



2(1 -a^) 25 - 272a + 25a^ + 12(l-8a + a2)lna ln(l - a) 



2a(12 - 18a - 236a^ + 19a^ - 48a7r^) In a - 12/i(a) In 



cm}, (B.16) 



+ 24(-l + 8a - 30a2 - 16a=^ + 2a^ - 12aMna) L2(a) + 8640^ [13(0) 
which again is explicitly UV and IR finite. 

B.7 Diagram XII 

For diagram XII only a four-particle cut is possible. Using the same techniques as before 
one finds: 



ImXII 



f 

Jm 



ds 



mlp^\s)-2pt>{s) 



(1), 



A{ml - sf 



X < In 



ml - s 1 , 



mt 



m 



- In -3 K «(2^6 + 5s) In ^ - 2(^6 - - 3^) 



2 2 m 



-2sml{ml + 3s) In ^ + !!^L^(s2 + 28sm2 - 5m^) 
s 3 
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2 fi 

184327r' 



;{l37 - 1024a - 324a^ + 1216a^ 



96a^ TT^ 



-24/i(a) 



21n(l 



a) — In — 



- 12a^(102 + 8a + 5a^) In a - 576a^ [In a ln(l - a) + L2(a)] }. (B.17) 

The IR divergent logarithm of A cancels the IR divergence in the sum of the diagrams II 
and 11^, thus rendering an IR finite result. Similarly, the dependence on the renormalization 
scale n cancels in the sum II + 11^" + V + XI + XI^ + XII, as it should be for the semileptonic 
decay rate T{b ^ utu). 

It is important to emphasize that these results do not change under Fierz-transforma- 
tions. This has been checked explicitly for the seven diagrams discussed in this appendix. 



ImX + XT =^^—l\ ^v]-^ + 4in^ 

L ^ J|m,=o 768 I 12 /i2 



5 



Im[XI + XI"^] 
In the notation of B this reads: 



^^/13_ 2_^^!!i\ 

mc=o 768 I 3 ^ ^ ' 



^12 /i2 



55 ml 



^ 139 2 ml 

Gd = -— + 7r2 + 41n^, 
12 jj,^ 

^ 13 2 1 

Gf = TT^ -In^. 

^ 3 fi^ 



B.8 The Limit ^ 

Here, we show that our results are well-defined in the limit and compare them 

with the corresponding quantities obtained in 0. From (|B.7|) , ( p.lO|) , ( [B.13|) and 
we get: 

iniViL^^o 768 1^12 /i^ j ' 

Im Villi o = ^f-^+41n!^V 
'™==° 768 V 3 fx^ J' 
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These results are scheme-dependent. In order to compare with [Q, where the 't Hooft- 
Veltman scheme is used, we need to add the corresponding matching coefficients, B^, B^, 
Bf,Bg evaluated using NDR. Again in the notation of one has: 

Bd = B, = 5, Bf = Bg = ^ (B.26) 

with B = Bd + Bf. + Bf + Bg = 11, cf. ( p.lO| ). The agreement with for the scheme- 
independent quantities Gi + Bi is evident, which provides a non-trivial check of our cal- 
culation. 



Appendix C 



As mentioned in the text, the calculation of Gb and Gd is much simplified by using the 
results for the spectral densities of the vector correlation function, l/vrlmll^,^ = q^qu pX + 
g^v pXi Eq. ( 2.22|) , which can be obtained from [|1^] when the running mass is replaced by 
the pole mass; we have also calculated the spectral densities directly by using Cutkosky 
rules. Defining 

pr = 3(pf)+c^pr +...), (c.i) 

where p^^^ is the tree-level and pf^ the next-to-leading order contribution, we find: 



pS\s) 



127r2 



s 



1 + 



2ml 



127r2 



(C.2) 
(C.3) 



9l 



1927r^ 



4 1- 



2ml 



Inl^ln fl-^~ 
s V s 



+ 2U 



m„ 



-4 1 + 2^-2^ In ^-4 1 ^ 1 + 5— 

\ 5 SIS \ S / \ S 



3mt \ , s — 



2s2 



In 



+ 3 1 



m^ 
1 + 3^ 
s 



16m^ 
3s2 



(C.4) 



2s 



5ml 
6^ 



(C.5) 



We repeat that the above spectral densities are expressed in terms of on-shell quark masses. 



In the hmit mc — > one obtains 



massless 



9l 



massless 



647r4 



(C.6) 



and lifj^u becomes transverse, as it should. 

Since for the calculation of diagram V, App. B.2., the part of pi corresponding to the 

(2) 

light quark self-energy is needed separately, we give in addition to the sum, pl , also the 
contributions of the single diagrams. All the three graphs are calculated in Feynman gauge 
and expressed in terms of the running MS mass whose relation to the pole mass is given 

by 



47r2 



(C.7) 



We denote the spectral density of the light quark's self-energy diagram by p^^^, that of 
the heavy quark by and the gluon-exchange diagram by pf^. 



11527r4 



Jn + 17!^-34 



— 6 In 



+ 6 



' s — ^ 2ml + s 



, "^r ^1 s — ml 

In ^ + 21n 

/^^ m^ 



(C.8) 
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23 047r4 



(s - m^) ^ 28 - 5 ^ - 17 f \+ 6{2s - 3ml) In ^ 



' s — ml\ ,^ 2 \ fi "^r oi s — ml 
^ {2ml + s) i In ^ + 2 In 



P^ 



(C.9) 



9l 



11527r4 



m„ 



11-55 



362 



30 In 



mt 



„s — ml I (f^lS^ 
+ 6 ^ < 1 + ^ + 34 ^ 



In 



p^ 



-24 



m,. 



m:. 



In 
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23 047r4 



[s - mi) ^ 28 + 157 ^ + 181 l"^) 



- 6(10s + 3ml) In 



mt 



vx 1 ( fn'i\'^\ , mi 

_6(.-,„^)^2 + 17^ + 17U In^ 



+ 24 



s — mt. 



(2s + m^) In 



s — m^ 
ml 



(C.ll) 



^2 s - 



2887r4 s 



'2 5 + 7^ 
s I \ s , 



In 



s — mt. 



mt 



-3 1 + 2^^ 



ln!!| + 21n4ln^-4L. f 



+ 12^ In 4 

s mi 



(C.12) 



9l 



11527r4 



-(s - m^) <( 46 - 23 ^ - 5 f ^ V 1 + fl + 2 ^ 1 In 



+ 6 



s — mt 



2{s + 2m^) In + (2s + 



mt. 



1 "^c 

^J? 



+ 21n4ln^-4L, f 

m2 s \ s ^ 



(C.13) 
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Figures 




Figure 1: The diagrams contributing to the forward-scattering amphtude Eq. (|2.16|) up to 
order as- The crossed circles denote insertions of the operator O. Of the three internal 
quark lines, the upper one denotes the c quark, the lower one the d quark, and the middle 
one the u antiquark. 
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Figure 2: Illustration of the effect of Fierz-transformations on the diagrams VI and XI. In 
contrast to the left-hand side diagrams, the ones on the right-hand side are well defined 
in NDR. 

Fig. 3 
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\sqrt{a} = mc/mb 



Figure 3: Radiative corrections from the different combinations of operator insertions as 
defined in Eq. ( |0|) as functions of mc/nih for /i = m^: short dashes: Cn, long dashes: 
Ci2, solid line: €22- The grey bar indicates the range of realistic values of mc/mh. The 
numerical values are tabulated in Tab. ^. 
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Fig. 4(a) 
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Fig. 4(b) 
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0.2 0.4 0.6 0.8 1 
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Figure 4: (a) The next-to-leading order correction /(a,/i) to Vih ceu), Eq. ( ^.41) , as 
function of i/a = mc/mh, normalized to one at = 0, for ^^(m^) = 0.117 and three 
different choices of the renormalization scale: solid line: fi = rrih, long dashes: /i = mb/2, 
short dashes: /i = 2mfe. (b) The next-to-leading order correction J{a,fj,) to T{b — > cud), 
Eq. ( |5.2| ), as function of ^/a = rric/mb, normalized to one at rric = and using the same 
parameters as in (a). The grey bar indicates the range of realistic values of rric/mb. 
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Figure 5: The semileptonic branching ratio of the B meson including non-perturbative 
corrections as a function of rrib with — rric fixed by HQET and as{mz) = 0.117. The 
three lines correspond to three different choices of the renormalization scale: solid line: 
/i = nib, long dashes: /i = mb/2, short dashes: /i = 2mb. 
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Tables 



Diagram No. 








I 


3 


1 


3 


II 


3Cf 


Cf 


3Cf 


III 


3Cf 


Cf 


3Cf 


IV 


3Cf 


Cf 


3Cf 


V 


3Cf 


Cf 


3Cf 


VI 


3Cf 


Cf 





VII 





Cf 


3Cf 


VIII 





Cf 





IX 





Cf 


3Cf 


X 





Cf 





XI 


3Cf 


Cf 





XII 


3Cf 


Cf 


3Cf 



Table 1: The colour-factors multiplying the diagrams shown in Fig. |I] for all possible 
combinations of operator insertions. For an arbitrary number of colours N^-, Cf equals 
(iVc'-l)/(2iVc). 








0.1 


0.2 


0.3 


0.4 ^ 0.5 


Ha{a) 


-3.62 


-3.22 


-2.71 


-2.21 


-1.75 




1.5 


2.43 


4.20 


6.70 


10.5 ^ + |7r2-151n2-31n(i-a) 




-3.62 


-3.58 


-3.44 


-3.11 


-2.16 


Hdia) 


1.5 


3.37 


6.46 


12.8 


24.3 



Table 2: Numerical values of the functions defined in Sec. 4. 
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y^ = mc/mb 0.1 0.2 0.3 0.4 

T{b^ccs) 1.112 1.005 0.671 0.283 0.043 

x6An^/{Glml) ±0 ±0.006 ±0.008 ±0.005 ±0.001 

Table 3: The decay rate r{b —>■ ccs) in units of G'j^ml / {GAn^) as function of mc/mb for 
/i = mi, and as{mz) = 0.117. The error denotes the uncertainty in the radiative corrections 
to the free quark decay. 



i/a = rric/ nib 


Cll 


Cuifi = nib) 


C22 









31 _ 


0.1 


-1.2 


-11. 


-1.8 


0.2 


0.32 


-11. 


-1.3 


0.3 


2.0 


-10. 


-1.0 


0.4 


3.8 


-9.5 


-0.73 


0.5 


5.6 


-8.9 


-0.51 


0.6 


7.6 


-8.4 


-0.33 


0.7 


9.9 


-7.8 


-0.20 


0.8 


13. 


-7.2 


-0.094 


0.9 


17. 


-6.6 


-0.027 


1 




-6 






Table 4: Numerical values of the coeffiencts qj defined in ( |5.1| ). 



as{niz) 


Parton Model 


nj HQE ig HQE [this work] 


0.110 


0.132 


0.130 0.121 


0.117 


0.128 


0.126 0.116 


0.124 


0.124 


0.121 0.111 



Table 5: B{B Xev) in different models depending on as{niz)- Input parameters: 
nib = 4.8 GeV, m^. = 1.3 GeV, which corresponds to Ai = —0.6 GeV^. In the phase-space 
factor of r(6 ccs) we use = 0.2 GeV. HQE is short-hand for heavy quark expansion. 
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